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We study the first observed charmless B — * VV modes, the B — > <j)K* decays, in perturbative 
QCD formalism. The obtained branching ratios B(B — ^ 4>K*) ~ 15 x 10 -6 are larger than ~ 9 x 10 -6 
from QCD factorization. The comparison of the predicted magnitudes and phases of the different 
helicity amplitudes, and branching ratios with experimental data can test the power counting rules, 
the evaluation of annihilation contributions, and the mechanism of dynamical penguin enhancement 
in perturbative QCD, respectively. 



The branching ratios of the penguin-dominated B — » Kir decays, about 3-4 times larger than those of the tree- 
dominated B — > 7T7T decays, indicate that penguin contributions must be enhanced. This enhancement can be achieved 
I— j | either by large Wilson coefficients C4.6 associated with the penguin operators in perturbative QCD (PQCD) or 
by a large chiral symmetry breaking scale mo associated with the kaon in QCD factorization (QCDF) The latter 

mechanism, called chiral enhancement, corresponds to a characteristic scale of 0(nib), at which we have mo(mt) ~ 3 
GeV and the smaller Wilson coefficients C4 e( m b)- The former mechanism, called dynamical enhancement, corresponds 

to a characteristic scale of 0(y Amj), A = Mb — rrif, being the B meson and b quark mass difference, at which we have 
C^ ' m (y/Ami,) ~ 1.5 GeV and the larger Wilson coefficients ^(v^Amj) ~ l.hCu^rrn,). Recently, we have proposed 
\^ \ the B — > 4>K decays as the appropriate modes to clarify the above issue ||f7| . These modes are not chirally enhanced, 
y—i • because <f> is a vector meson, and insensitive to the variation of the unitarity angle 03 , because they are pure penguin 

processes. If the data of the branching ratios B(B — > <f)K) are settled down at values around 10 x 10~ 6 instead 

of 4 x 10 -6 |s|,|lO|], the dynamical enhancement of penguin contributions to charmless nonleptonic B meson decays 

will gain a strong support. 

Here we argue why the characteristic scale involved in two-body B meson decays must be of 0(\/AMb) in PQCD 
' from two points of view. Consider a two-body nonleptonic decay, in which the two final-state light mesons move 
back-to-back with large momenta. The lowest-order diagram for its amplitude contains a hard gluon attaching the 
O ' spectator quark. Intuitively, the spectator quark in the B meson, forming a soft cloud around the heavy b quark, 
carries momentum of order A. The spectator quark on the light-meson side carries momentum of O(Mb) in order to 
1-^ ' form the fast-moving light meson with the u quark produced in the b quark decay. Note that the end-point singularities 
from the small spectator momentum on the light-meson side do not exist in a self-consistent PQCD formalism, because 
of Sudakov suppression from and threshold resummations Based on the above argument, the hard gluon is 

off-shell by order of AMb- This scale characterizes the corresponding quark- level hard amplitude, which involves the 
four-fermion decay vertex. Theoretically, the hard scale AMb is essential for constructing a gauge invariant B meson 
wave function. This wave function, though being a nonlocal matrix element, is gauge invariant in the presence of 
the path-ordered Wilson line integral. A careful investigation Jl3|,[l4| shows that the 0(a 2 s ) diagram with the second 
gluon attaching the hard gluon contributes to this line integral. That is, this diagram contains the soft divergence, 
which is factorized into the B meson wave function. This is possible, only when the hard gluon is off-shell by the 
intermediate scale AMb, rather than by M%. 

In this work we shall perform a PQCD analysis of the first observed charmless B — > VV modes, the B — > cj>K* 
decays, which are, similar to B — > (f>K, also appropriate for distinguishing the different penguin enhancing mechanism. 
Besides, the B — > VV modes reveal dynamics of exclusive B meson decays more than the B — > PP and VP modes 
through the measurement of the magnitudes and the phases of various helicity amplitudes. According to the power 
counting rules defined in |?]] , the longitudinal amplitude is leading, and the other two amplitudes are down by a power 
of M^/Mb or of Mk*/Mb, M<f, and Mk* being the cf> and K* meson masses, respectively. Since the B — * <f>K* 
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decays are insensitive to the unitarity angle, the relative phases among the helicity amplitudes mainly arise from 
strong interaction. The annihilation contributions, which can be evaluated unambiguously in our approach, generate 
the strong phases. Therefore, comparing the predicted magnitudes and relative phases among the different helicity 
amplitudes, and the predicted branching ratios with experimental data, we test the power counting rules, the evaluation 
of annihilation contributions, and the mechanism of dynamical penguin enhancement in PQCD, respectively. 

The idea of the PQCD factorization theorem for two-body nonleptonic B meson decays has been reviewed in 
p|, p^|Jl^ ], which is subject to corrections of 0{a 2 ) and 0{K/Mb)- In this formalism decay amplitudes are expressed 
as the convolutions of the corresponding hard parts with universal meson distribution amplitudes [l^Jl^| , which are 
regarded as the nonperturbative inputs. Because of the Sudakov effects from kr and threshold resummations, the 
end-point singularities do not exist as stated above. Therefore, PQCD involves inputs less than in QCDF, for which 
form factors, meson distribution amplitudes, and infrared cutoffs for regulating the end-point singularities are all 
independent parameters Strictly speaking, the infrared cutoffs, signifying important soft contributions to the 

nonfactorizable and annihilation amplitudes, imply that the factorization formulas in QCDF are not self-consistent. 

We work in the frame with the B meson at rest, i.e., with the B meson momentum P\ — (M B /V2)(1, 1> Or) in the 
light-cone coordinates. Assume that the <p (K*) meson moves in the plus (minus) z direction carrying the momentum 
Pi (P3) and the polarization vectors e 2 (€3). The B — > (j)K* decay rates are written as 

n _ G\P C ^ M W M W ; (1) 



16ttM| 



B „- 



L,T 



where P c = I-P22I = \Paz\ 1S the momentum of either of the outgoing vector mesons, and the superscript a denotes 
the helicity states of the two vector mesons with L(T) standing for the longitudinal (transverse) component. The 
amplitude M.^ is decomposed into 



a cT + P?Pi + i ^ vaS) Pi a Pia 

y M^Mk* 1 1 M^M K * . 



= M 2 B M L + M 2 B M N e*(a = T) ■ e* 3 (a = T) + iM T ^ lp ^ a {oY^{a)P^P 3p , (2) 



with the convention e 0123 = lQ and the definitions, 



Ml M L = a e* 2 (L) ■ e* s (L) + \ e* 2 (L) ■ P 1 e* 3 (L) ■ P 1 , 

M% M N = a e* 2 {T) ■ e%{T) , (3) 
\l, 



We define the helicity amplitudes, 



M<j,M K * 



A a = -iM 2 B M L , 
Ay = tV2M B M N , 



A± = {M^Mk. ^2{r 2 - l)M T , (4) 



with the normalization factor £ = \J Gj^Pc/ (16ttM b T) and the ratio r = P 2 -P 3 j (M^Mk* ) ■ These helicity amplitudes 
satisfy the relation, 

L4 | 2 + |A||| 2 + L4J 2 = 1, (5) 
following the helicity summation in Eq. (|l|). We also introduce another equivalent set of helicity amplitudes, 

H Q = M 2 M L , 



H ± = M 2 M N T M 4 M k , ^r 2 - IM T , (6) 

with the helicity summation, 



1 This convention corresponds to tr(js fi fi fc, fl) = —4ie al3 ~ /p a a bpc~ f d p . 



2 



Y^M^M^ = \H \ 2 + \H+\ 2 + \H-\ 



(7) 



The B — > <j)K* decays involve the emission and annihilation topologies, both of which are classified into factorizable 
diagrams, where hard gluons attach the valence quarks in the same meson, and nonfactorizable diagrams, where hard 
gluons attach the valence quarks in different mesons. The amplitudes are written as 



m h = Uv?f£ + v;m%\ + f B V t *Ff a + V t *Mf a , 



M H = 



V t *M%l + f B V?F$ + V t *M%l - f B V:F Ha - v:M Ha , 



(8) 
(9) 



for the B° 



4>K*° and B + — ► <fiK* + modes, respectively, where the subscript H = L,N,T denotes the different 
helicity amplitudes, e (a) denotes the emission (annihilation) topology, and V q — V* s V q b are the products of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The hard parts for the factorizable amplitudes F and for 
the nonfactorizable amplitudes M. are derived by contracting the following structures to the lowest-order one-gluon- 
exchange diagrams, 



= (/\ + ^b)75$M) , 



M, 



[M K , ML)$k-(x)+ ML) fo& K .{x) + M K ,I<S> s K ,{x)] , 

M K 



K y „p*'K'f%{T)P$n°_*%{x) 



M K . MTWk*(x)+ MT) ftQie-ix) 



Pa-n, 



ie^ p<7l5 r^(T)P 3 p n^ a K ,(x) 



(10) 

(11) 
(12) 

(13) 
(14) 



where n + = (1,0, Ot) and n_ = (0, 1,0t) are dimensionless vectors on the light cone. Equations ( |l l| ) and ( |12| ) are 
associated with the longitudinally and transversely polarized <f> mesons, respectively. The structures associated with 
the K* meson are similar as shown above. 

To extract the contributions to the helicity amplitude Ml, the following parametrization for the longitudinal 
polarization vectors is useful: 



ML) = 



El. 

M A 



Ms 



Pi-n- 



M L ) = 



M, 



K* 



M, 



K* 



P* ■ n. 



-n+ 



(15) 



which satisfy the normalization e\{L) = e|(I>) = — 1 and the orthogonality €2(L) ■ P^ = (j,{L) • P3 = for the on-shell 
conditions P% = M? and Pf = Mjj-,. We first keep the full dependence on the light meson masses M$ and Mk* in 
the momenta P2 and P3. After deriving the factorization formulas, which are well-defined in the limit M^Mk* — > 0, 
we drop the terms proportional to r?, r 2 K , ~ 0.04, with the ratios — M^/Mb and tk- = Mk*/Mb- Under this 
approximation, the expressions of the 



and K* meson momenta are then as simple as 



Po 



M B 



(1,0, Or), Ps 



M f 



(0,1, Or) 



(16) 



For the extraction of the helicity amplitudes Mn and Mr, Eq. JL6| ) and the transverse polarization vectors, 

e 2 (T) = (0,0,l T ), e 3 (T) = (0,0,l T ), (17) 

can be adopted directly. The explicit factorization formulas are collected in the Appendix. 

The power counting rules in PQCD |Q tells that the factorizable amplitude Fl c (corresponding to the B — > K* 
transition form factor) is leading, and the other factorizable amplitudes are at least down by a power of or r^. . 
The nonfactorizable amplitudes M are suppressed by a power of A/ Mb- Hence, the formalism presented in this work 
is complete at 0(M^^k* /Mb), and subject to corrections of 0(A/Mb)- Equation (||) then implies that the helicity 
amplitude Aq is leading in the heavy-quark limit, and An and A± are next-to-leading. The factorizable annihilation 
amplitudes Fn a , being suppressed only by M^^k* /Mb and almost imaginary, are the major source of the strong 
phases in PQCD. Since the B — > <f>K* decays are the pure penguin processes with a weak dependence on the unitarity 
angle ^3, these strong phases determine the relative phases among the helicity amplitudes Aq, An and A±. 
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For the B meson wave function, we employ the model [|l| , 



<5>b(x, b) = N B x 2 (l - xf exp 



1 (xM B \ 2 u%b 2 



2 V LO B J 



(18) 



where the shape parameter ljb = 0.4 GeV has been adopted in all our previous analyses of exclusive B meson decays. 
The normalization constant Nb = 91.784 GeV is related to the decay constant fs — 190 MeV (in the convention 
Jtc = 130 MeV). It is known that there are two B meson wave functions and <f>B, which are related to the three- 
parton B meson wave functions through a set of equations of motion |L7] ^0) . Because of the unknown three-parton 
wave functions, the equations of motion in fact do not impose any constraint on the functional form of $s and $b- 



Our simple choice of the model wave functions corresponds to <f>s in Eq. ( |18| ) and $s — 0. This choice is legitimate, 
since the contribul 
current formalism. 



since the contribution from <&b is suppressed by a power of A/ Mb M], and negligible within the accuracy of the 
The </> and K* meson distribution amplitudes up to twist 3 are given by pl[ 
3/0 



<& K ,{x) 
& K .(x) 

& K .{x) 



f2N c 
f T 



2V2N, 

f T 

J<t> 



x(l — x) , 

3(1 - 2a;) 2 + 1.68C] /2 (1 - 2x) + 0.69 



1 + (1 - 2x) In 



1-x 



4V2V, 
3/J 

y/2W c 

S<t> 



3(1 - 2s) (4.5 - 11.2a; + 11.2a: 2 ) + 1.38 In 



x{l-x) l + 0.2C2 /2 (l-2a: 
3 



1-x 



2^J2N r I 4 
3/ 



~[1 + (1 - 2a;) 2 ] + 0.24[3(1 - 2.x) 2 - 1] + 0.96C 4 1/2 (1 - 2x) 



<{^/2N l 

\f2N, 
Jk- 



(1 -2a;) 1 + 0.93(10a; 2 - 10a; + 1) 
x{l -x)\l + 0.57(1 - 2a;) + 0.07C2 /2 (1 - 2x) 



2J2N r 



0.3(1 - 2x) [3(1 - 2xf + 10(1 - 2x) - lj + 1.68C 4 1/2 (1 - 2x) 



0.06(1 - 2a;) 2 [5(1 - 2a;) 2 - 3] + 0.36 [1 - 2(1 - 2a:) (1 + ln(l - x))} j 
3(1 - 2x) [l + 0.2(1 - 2x) + 0.6(10a; 2 - 10a; + 1)1 



2y/2N, 

-0.12x(l - x) + 0.36[1 - 6a; - 21n(l - x)} 
3/£. 



IK* 

/2N, 
Ik* 



-.x{\ -x) 1 + 0.6(1 - 2a:) + 0.04C| /2 (1 - 2a; 



1 + (1 - 2x) 2 + 0.44(1 - 2a;) 3 + 0.4C 2 1/2 (1 - 2x 



2^J2N r \^ 
+0.88C* 4 1/2 (1 - 2x) + 0.48[2a; + ln(l - x)\ 
Ik 



iV2N c 



|3(1 -2x) 1 + 0.19(1 - 2x) + 0.81(10a; 2 - 10a; + 1) 



-1.14x(l -x) + 0.48[1 - 6a; - 21n(l - a;)]} 
with the Gegenbauer polynomials, 



Cl /2 {0 = \('Se - 1) , C 4 1/2 (O = i(35£ 4 -30£ 2 + 3), C 2 3/2 (0 = ^(5C 2 -1) 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 

(26) 

(27) 
(28) 

(29) 
(30) 
(31) 
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We employ G F = 1.16639 x 1CT 5 GeV -2 , the Wolfenstein parameters A = 0.2196, A = 0.819, and R b = 0.38, the 
unitarity angle 03 = 90°, the masses Mb — 5.28 GeV, = 1.02 GeV and Mk- = 0.89 GeV, the decay constants 
/ = 237 MeV, fj = 220 MeV, f K * = 200 MeV, and /£, = 160 MeV, and the B° d (B+) meson lifetime t b o = 1.55 
ps (t b + — 1.65 ps) j22|. We have confirmed that the above distribution amplitudes and decay constants lead to the 
B — > K* transition form factors |23] in agreement with those from light-cone Q CD s um rules |24j. We have also 



confirmed that the averaged values of the running hard scales t defined by Eqs. ( A2C ) and ( |A2lj ) in the Appendix 



are indeed about \J AMg ~ 1.6 GeV. Note that the B — > <j>K* branching ratios are insensitive to the variation of 
<f>3. The results for the helicity amplitudes Aq, A\\ and A±, including their relative phases </>|| = Arg(A\\/Ao) and 
4>±_ = Arg(A±/Ao), are displayed in Table I. The contributions to the B — > 4>K* branching ratios mainly arise 
from the longitudinal polarizations Aq, because of the relation \Ao\ 2 ^> \A^\ 2 ~ |^4^| 2 , which is expected from 
the power counting rules. It is easy to observe that the ratios \H-/Hq\ 2 and \H+/Hq\ 2 obtained in PQCD are 



close to those in QCDF 25 . The annihilation contributions are the major source of the strong phases, and the 
nonfactorizable contributions are the minor one. The values of 0|| and (f>± in the rows (I)- (III) of Table II indicate 
that the phases from the former are about 4-5 times those from the latter (but opposite in sign). Without these 
sources, we have <f>n = <fi± = w. Note that the relative phases among the different helicity amplitudes can not be 
predicted unambiguously in QCDF, due to the arbitrary complex cutoffs for the evaluation of the nonfactorizable and 
annihilation contributions. 

We examine the theoretical uncertainty from the variation of the hard scales t, which are defined as the invariant 
masses of the internal particles, and required to be higher than the factorization scales 1/6, b being the transverse 
extents of the mesons. This examination estimates higher-order corrections to the hard amplitudes, which are the most 
important theoretical uncertainty for penguin-dominated B meson decays. The light meson distribution amplitudes 
have been determined in QCD sum rules. The possible 30% variation of the coefficients of the Gegenbauer polynomials 
in these distribution amplitudes lead only to little changes of our predictions. We consider the hard scales t located 
between 0.75-1.25 times the invariant masses of the internal particles. The predictions for the B — ► <f>K branching 
ratios from the above range are consistent with the data with uncertainty Q]. We then obtain the B — > 4>K* branching 
ratios, 

B(B° d -» 4>K*°) = (14.86^-88) x 10~ 6 , B{B ± -> QK**) = (15.96t|;£) x 10~ 6 . (32) 

The relative phases 4>\\ and (j>±, and the magnitudes |^4o| 2 j l^||| 2 an d |^4._l | 2 of the helicity amplitudes are quite stable 
under the variation of the hard scales t. They change within 0.05 rad. and within 0.01, respectively. There is another 

(T) (T) 

minor source of theoretical uncertainty from the light meson decay constants and f K J- If they reduce by 5%, the 
predicted branching ratios will decrease by 10%. The CP asymmetries of the B — * (j>K* modes are, as of B — > <f>K, 
vanishingly small (less than 2%). 

The above branching ratios are larger than those from QCDF p5[, 

B(B° d -> <j)K*°) = 8.71 x 10~ 6 , B(B ± -> c/jK*^) = 9.30 x 10~ 6 , (33) 

due to the dynamical enhancement of penguin contributions. We emphasize that the annihilation amplitudes, though 
not negligible, are not responsible for the large branching ratios in PQCD, since they are mainly imaginary. This 
is understood by comparing the branching ratios in Table I and in the row (II) of Table II. The nonfactorizable 
contributions are not either as shown by the branching ratios in Table I and in the row (III) of Table II. However, the 
annihilation contributions, parametrized as being real, are important in QCDF in order to explain the large B — > <pK 
branching ratios. With the almost real annihilation contributions, the B — > tpK branching ratios obtained in QCDF 
can increase from 4 x 10~ 6 to 7 x 10~ 6 M. The values quoted in Eq. ( |33| ) do not include the annihilation contributions. 
The current experimental data of B(L P — > cf>K* 



CLEO [26] 
BELLE H] 
BABAR H 

and those of B{B ± -> (pK^), 

CLEO [26] 
BELLE §7| 
BABAR H 



(11.5±| : f±i;f) x 10- 6 , 
(15±| =b 3) x 10~ 6 , 

(8.6±H ± L1 ) x 10_6 > ( 34 ) 



Cio.6±£S±i5) >< io- 6 , 

< 36 x 10~ 6 , 

(9.7if:4 ± 1.7) x 10 -8 , (35) 
are not yet precise enough to distinguish the two different approaches. 
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In this paper we have studied the first observed B — > VV modes, the B — > </>if* decays, using the PQCD formalism. 
It has been stressed that two-body heavy meson decays are characterized by a scale of 0(AMb) in PQCD, for which 
penguin contributions are dynamically enhanced. This enhancement renders penguin-dominated decay modes acquire 
branching ratios larger than those in QCDF, even when the final-state particles are vector mesons. We have proposed 
the B — > (f)K^*' decays as the ideal modes to test the significance of this mechanism. If their branching ratios are as 
large as 10 x 10~ 6 (15 x 1CP 6 ) (independent of the unitarity angle ^3), dynamical enhancement will be convincing. 
We have also emphasized that the relative importance and the relative strong phases among the different hclicity 
amplitudes in the B — > VV modes can be predicted unambiguously in PQCD, which are determined by the power 
counting rules and by the annihilation contributions, respectively. These predictions are insensitive to the variation 
of the hard scales. Therefore, the comparison of the results presented here with future experimental data will provide 
a stringent confrontation of the PQCD approach. 
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APPENDIX A: FACTORIZATION FORMULAS 



In this appendix we present the explicit expressions of the factorizable and nonfactorizable amplitudes in Eq. 
The effective Hamiltonian for the flavor-changing b — > s transition is given by 



i=3 



with the CKM matrix elements V q = V* s V q b and the operators 

0[ q) = (siqj)v-A(qjbi)v-A , = (siqi) V -A(qjbj) V -A , 

3 = {sibi)v-A^(q q])v-A , 4 = {sibj) v ^A^2{q]qi)v-A , 
9 9 

1 q 

°7 = -(sibi)v-A22e q (qjqj) v+A , 8 = -(sibj)v-A ^2 e i%n)v+A 
9 9 

3 



O9 = -(sibi) v ~A ^2 e q(lj1j)v-A , O w = -(sibj)v-A ^ e <3 (g i q i ) 



V-A 



(Al) 



(A2) 



i and j being the color indices. Using the unitarity condition, the CKM matrix elements for the penguin operators 
O3-O10 can also be expressed as V u + V c = —Vt. The unitarity angle <p3 is denned via 



Kb = |Kfe|exp(— i<j> 3 ) . 
Here we adopt the Wolfenstein parametrization for the CKM matrix upto 0(A 3 ), 



V u d V us V u b ' 

v cd v cs v cb 

. Vtd Vt a V t b , 



with the parameters |2t 



l-A 2 /2 A A\ 3 (p-ir])> 

-A 1 - A 2 /2 AX 2 

,AX 3 (l- p-irj) -A\ 2 1 



A = 0.2196 ±0.0023 
A = 0.819 ±0.035 , 



(A3) 



(A4) 



R b = y/p 2 + rj 2 = 0.41 ± 0.07 



(A5) 



G 



The factorizable amplitudes Fg } e and Fg } a = F^ 4 + are written as 

pi pOO 

F^=8ttC f MI d Xl dx 3 feid6i6 3 d6 3 $ B (xi,6i) 
Jo Jo 

x { [(1 + x 3 )$ K * (x 3 ) + r K , (1 - 2x 3 )(& K , (x 3 ) + $ S K , (x 3 ))] 

+2r K ^ s K ,(x 3 )E^(t^)h e (x 3 ,x u b 3 ,b 1 )} , (A6) 

pi pOO 

F { N q) e = 8ttC f M 2 b d Xl dx 3 hdhhdh^Bixuh) 
Jo Jo 

xr^fCfe) + 2r K *$ v K ,(x 3 ) + r K *x 3 ($ v K ,(x 3 ) - $ a K .{x 3 ))] 
xEi"\t^)h e {x u x 3 MM) 

+r K ,[^ v K ,(x 3 ) + ^ a K ,(x 3 )}E^(t^)h e (x s ,x 1 ,h,b 1 )} , (A7) 

pi poo 

F$=16nC F Ml d Xl dx 3 bidbib 3 db 3 <$> B {xiM) 
Jo Jo 

xr4K,(x 3 ) + 2r K ,<P a K ,(x 3 ) - r K ,x 3 (^ K ,{x 3 ) - $ a K ,{x 3 ))] 
xE^(t^)h e (x u x 3 ,b u b 3 ) 

+r K 4^ K 4x 3 ) + ^> a K ,(x 3 )}E^(t^)h e (x 3 ,x u b 3 ,b 1 )} , (A8) 



pi pOO 

F { l ^ 4 ^8ttC f MI / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x { [-(1 - x 3 )<5>^x 2 )<5> K * (x 3 ) + 2r 4> r K *$ s 4> (x 2 ){x 3 $ t K , (x 3 ) + (2 - x 3 )$^, (x 3 ))] 
xE { £{t^)h a {x 2 ,l-x 3 MM) 

+ [x 2 $ (x 2 )$ K .(x 3 ) + 2 r<p r K ,^ K ,{x 3 )((l - x 2 )^%{x 2 ) - (1 +x 2 )%{x 2 ))] 

x E^(t^)h a (l-x 3 ,x 2 ,b 3 ,b 2 )} , (A9) 

pi POO 

F nL = -%kC F M 2 b / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x r<p r K , { [(2 - x 3 ){*l(x 2 )& K . (x 3 ) + ^(x 2 )$ a K , (x 3 )) 
+x 3 (*l(x 2 )* a K . (x 3 ) + %{x 2 W K , (x 3 ))] E^(t^)h a (x 2 , 1 - x 3 , b 2 , 63) 
- [(1 + z 2 )(<I>£(x 2 )<I>^ (a*) + *5(a:2)*S-. (* 3 )) 

-(l-^ 2 )($^(x 2 )$^(x 3 ) + $5(x 2 )<i>^(x 3 ))]^ 4 ) (4 2 ))/ la (l-x3,x 2 ,6 3 ,6 2 )} , (A10) 

= -167rC F M| / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x r^r K , { [x 3 (% (x 2 ) $ V K , (x 3 ) + % (x 2 ) $ a K , (x 3 ) ) 

+ (2 - x 3 )($l(x 2 )$ a K , (x 3 ) + $%{x 2 )& K ,{x 3 ))] E^(t^)h a (x 2 , 1 - x 3 ,6 2 ,6 3 ) 
+ [(1 - x 2 )(^(x 2 )^(x 3 ) + $%{x 2 )3> a K *{x 3 )) 

-(1 + x 2 )(<i>;(x 2 )$^ (x 3 ) + *%(x 2 )W K . (x 3 ))] E^(t^)h a (l - x 3 , x 2 , 63, b 2 )} , (All) 
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pi pOO 

F^g = lQirC F M% / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x { [r K , (1 - x 3 )$4x 2 )($* f , (a*) + $ S K , (x 3 )) - 2r 4> %(x 2 )$ K * (x 3 )] 
xE^(t^)h a (x 2 ,l-x 3 ,b 2 ,b 3 ) 

+ [r^x 2 (^(x 2 ) - *%(x 2 ))*k- (a*) + 2r K , ^(x 2 )^ s K , (x 3 )] 
xE { ^(t^)h a (l-x 3 ,x 2 ,b 3 ,b 2 )} , 

pi pOO 

f nL = IGnCpMl / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x {r*($J(s 2 ) + %{x 2 ))<S> T K , (x 3 )E^(t^)h a (x 2 , 1 - x 3 ,b 2 ,b 3 ) 
+ r K ^l(x 2 )(^ K ,(x 3 ) - $$(x 3 ))E$(tP)h a (l -x 3 ,x 2 ,hM} , 

pi pOO 

F^g = 32irC F M B / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x {r^l(x 2 ) + *2(s 2 ))*£. (X3)^ 6 ) (ii 1) )^(^, 1 - a*, 6 2 , 6 3 ) 

+ ^$J(x 2 )(<i>^(x 3 )-n*(^))4 9 6 ) (il 2) )^(l- ^3^2,6 3 ,fe 2 )} , 



(A12) 



(A13) 



(A14) 



The expression of the factorizable amplitudes Fn a from the tree operators Oi and 2 are the same as F^ aA but with 

the evolution factor E$ replaced by E$ . 

The factors E(t) contain the evolution from the W boson mass to the hard scales t in the Wilson coefficients a(t), 
and from t to the factorization scale 1/6 in the Sudakov factors S(t): 

E^(t) =a s (t)a^(t)S B (t)S K , (t) , 
Ej*> (t) =a s {t)a!f\t)S^t)S K ,{t). 
The Wilson coefficients a in the above formulas are given by 



(A15) 



a (q) 
"1 



,(?) 



,(9) 



C2 + aT 



3 ( n . C10 



C 6 



(1) , „(<?) 

tic 



C 7 



a 



&t resummation of large logarithmic corrections to the B, <j> and K* meson distribution amplitudes lead to the 
exponentials Sb, Sj, and Sk*, respectively, 



Sb (t) = exp 
S^t) = exp 
S K * (*) = exp 



- 2 



-s(xiP+,6i) 
-s(a; 2 P+,6 2 )-s((l 



1/61 M 



-x 2 )P+,b 2 )-2 f ^ 7 K(A 2 )) 



- S (x 3 P 3 -,63)- S ((l-X3)P 3 -,6 3 )-2 f ^ 7 ( as (/2 2 )) 

Ji/63 ^ 



(A16) 
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with the quark anomalous dimension 7 = —a s /iT. The variables 61, 62, and 63, conjugate to the parton transverse 
momenta k\x, k 2 x, and k 3 T, represent the transverse extents of the B, <f> and K* mesons, respectively. The expression 
for the exponent s is referred to |M)|]3l| ]. The above Sudakov exponentials decrease fast in the large 6 region [ ]lT| , |l2"|] , 
such that the B — ► <fiK* hard amplitudes remain sufficiently perturbative in the end-point region. 
The hard functions h's are 



h e (xi,x 3 ,h,b 3 ) 



h a (x 2 ,x 3 ,b 2 ,b 3 ) 



K (^/xix 3 M B bi) S t (x 3 ) 

x [6{h - b 3 )K {^M B bx) I (V^M B b 3 ) 

+9(b 3 - h)K (^M B b 3 ) I (^M B h)} , 

■ s 2 
ITT 

(i) 



(A17) 



(^xJM B b 2 ) S t (x 3 ) 



x 6{b 2 - b 3 )H^> {^M B b 2 ) Jo (V^M B b 3 ) 
+8(b 3 - b 2 )H^ (^M B b 3 ) J Q (y/x^M B b 2 ) . 

We have proposed the parametrization for the evolution function St (x) from threshold resummation |Tl] , |32| , 

2 1+2c r(3/2 + c), 



St(x) = 



ViT(l + c) 



[x(l-x)]< 



(A18) 



(A19) 



where the parameter c is chosen as c = 0.4 for the B — > (j)K* decays. This factor modifies the end-point behavior of the 
meson distribution amplitudes, rendering them vanish faster at x — > 0. Threshold resummation for nonfactorizable 
diagrams is weaker and negligible. Kq,Iq,Hq and Jo are the Bessel functions. 

The hard scales t are chosen as the maxima of the virtualities of the internal particles involved in the hard amplitudes, 
including 1/bf. 



4 1 ) = maxC^Ms, 1/6 X , l/6 3 ) , 

4 2 ) = max( y/xTM B , 1/h, l/b 3 ) , 

= max(VT — ~ 3 M B , l/b 2 , 1/63) , 
t< 2 > =max(V^M B , 1/63,1/63) • 



(A20) 



(A21) 



When the PQCD formalism is extended to 0(a 2 s ), the hard scales can be determined more precisely and the scale 
independence of our predictions will be improved. Before this calculation is carried out, we consider the variation of, 
for example, t e in the following range, 



max(0.75^3"M B , 1/6!, l/6 3 ) < < max(1.25^3"M B , 1/6!, l/6 3 ) , 
max(0.75VxTM s , 1/h, l/6 3 ) < tf> < max(1.25^/xlM B , 1/h, l/6 3 ) , 

in order to estimate the 0(ot 2 s ) corrections. The range for t a is chosen in a similar way. 



(A22) 



The nonfactorizable amplitudes Mxj e = -A^jjL 
on kinematic variables of all the three mesons IB,' 



] , are written as 



M3 L, depending 



_____ r^- r°° 

= lQirC F Ml^/W c I d[x] hdhb 2 db 2 $ B (xx,h) 
Jo Jo 

x {$<l>(x 2 ) [~{x 2 + x 3 )<l> K * (x 3 ) + r K *X 3 {& K , (13) + $° K . (x 3 ))] 



xEif(t^)h^( Xl ,x 2 ,x 3 ,h,b 2 ) 
+<P ct ,(x 2 ) [(1 - x 2 )$ K *{x 3 ) + r K *x 3 ($K*(x 3 ) - $Jr.(s a ))] 
xE { £\t { p)hf{ Xl ,x 2 ,x 3 ,h,b 2 )) , 



(A23) 



M 



(9) 



Ne3 - 16itC f M b ^2N c I d[x\ 

Jo 



hdhb 2 db 2 <& B (xi,bi) 



x r< t,{[x 2 (^(x 2 ) + ^(x 2 ))^(x 3 ) 
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-2r*. (x 2 + x 3 ){%(x 2 )* v K . (x 3 ) + %{x 2 )* a K . (xs))] 
xEif(t$ ) )hW(x 1 ,x 2 ,xs,h,b 2 ) 

+ (l-X 2 )(^(x 2 ) + %(x 2 ))^,(xs) 
xE^'(tf)hf(x u X 2 ,Xs,h,b 2 )} , 

pi pOO 

M^ ] e3 = 32 7 rC F Mlv / 2X / d[a;] / &id&i&2d&2*u(&i,&i) 

Jo Jo 



xr*{[z 2 ($2(a; 2 ) + *3(z 2 ))$£.(z3) 

-2r*. + ^x^)^. (x 3 ) + $;^ 2 )^. (x 3 ))] 

x^f (tf)h%\x u x 2 ,x 3 MM) 
+(l-x 2 )($l(x 2 ) + $>%(x 2 ))$ T K *(x 3 ) 

xE^'(tf)hf(x u x 2 ,x 3 MM)} , 

/' 1 pOO 

M { ZIs = 1&kC f M%^2N c d[x] b 1 db 1 b 2 db 2 <S>B(xi,b 1 ) 

Jo Jo 

xr { [-x 2 (<^0 2 ) - %(x 2 ))$K' (x 3 ) 
+r K 'X 2 (^(x 2 ) - %(x 2 ))(^ K , (x 3 ) - *' K .(x 3 )) 
+r K ,x 3 (^(x 2 ) + %(x 2 ))(& K . (x 3 ) + <f> s K , (S3))] 

x^'^^^i.ara.^.fti.fe) 

+ [-(1 - x 2 )(^(x 2 ) + ^(aa))**- (x 3 ) 

+TK* (1 - ^)(^(^2) + %(x 2 ))(& K . (x 3 ) - & K , (x 3 )) 

+r K ,x 3 ($%(x 2 ) - %(x 2 ))(^ K , (x 3 ) + $ s K .(x 3 ))] 
xEif(t^ ) )h^(x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , 

_ pi pOO 

J^Nes = -^C F M 2 B ^/W C / d[x] b 1 db 1 b 2 db 2 ^ B (xiM) 

Jo Jo 

xr K *xs$l(x 2 )($ v K ,(xs) - $ a K .(xs)) 

^{E ( £\tf)h { p(x 1 ,x 2 ,x 3 MM)+E ( £\tf)hf(x 1 ,x 2 ,x 3 MM)} 



M {q) - 2M {q) 



rl poo 

M^1 6 = -1&ttCfM 2 b ^/W c d[x] b 1 db 1 b 2 db 2 ^ B (x u b 1 ) 

Jo Jo 

x$0(x 2 ) {[x 2 <$> K ,(x 3 ) +r K ,x 3 (<P t K ,(x 3 ) - ^ s K .(x 3 ))] 
xElf'(t$ ) )hV(x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

+ [-(l-x 2 + x 3 )<P K , (x 3 ) + r K ,x 3 (¥ K , (x 3 ) + (&s))] 
xE < £'(tf)hf(x 1 ,x 2 ,x 3 MM)} , 

M { $ e(i = -IQKC F M 2 B ^/2N~ C d[x] b 1 db 1 b 2 db 2 ^ B (x 1 ,b 1 ) 

Jo Jo 
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xr^[x 2 (* J(x 2 ) - $$(z 2 ))*£. {x 3 )E { £' (4 1} )4 1} (x! , x 2 , x 3 , h, b 2 ) 

+ [(l-x 2 )($l(x 2 )-Z>l(x 2 ))$ T K ,(x 3 ) 

-2r K ,(l - x 2 + x s )($l(x 2 )^ K ,(x 3 ) - ^(x 2 )$ a K ,(x 3 ))] 

xE ( e f' (t^ ) )h^(x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , (A30) 

_ pi pOC 

M { j!l 6 = -327rC F M|V^Vc / / M&iM&2$b(zi,&i) 

Jo Jo 

xrt{x 2 ($l(x 2 )-^(x 2 ))*l.(x 3 )Eif' {t { p)h { p{ Xll x 2 ,x 3 MM) 

+ [(l-x 2 )(%(x 2 )-^(x 2 ))$ T K .(x 3 ) 

-2r K . (1 - x 2 + x 3 )(<t>l(x 2 )<P a K , (x 3 ) - *%{x 2 )& K . (x 3 ))] 

x^f(4 2) )^ 2) (^i^2,x 3 ,6 1 ,6 2 )} , (A31) 

_ P 1 /• OO 

M^ ) aS = W-KC F M 2 B ^/2N' c I d[x] hdhfodb^Bixuh) 

Jo Jo 

x {[(1 - a; 3 )$^(a; 2 )$if.(a;3) +r<t,r K * ((1 + x 2 - x 3 ){$> t 4> (x 2 )<$> t K , (x 3 ) 
-%(x 2 )^ s K ,(x 3 )) -(1-X2- x 3 )($l(x 2 )$ s K *(x 3 ) - *%(x2)* t K .(x 3 )))] 
xEi g 3 ) '(tf ) )hf\x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

- [x 2 ^4>( x 2)^k* (x 3 ) - 2r r K . ($* (a; 2 )^, (x 3 ) + %(x 2 )<P s K , (x 3 )) 
+r^r K , (1+X2- x 3 )(&+(x 2 )& K . (x 3 ) - %(x 2 )$ s K , (x 3 )) 
+ r<p r K * (1-X2- x 3 )(<S>^(x 2 )<f> s K , (x 3 ) - %{x 2 )& K , (x 3 ))] 

*E^\tf)hf (x 1 ,x 2 ,x 3 ,b 1 , b 2 )} , (A32) 

pi pOO 

M ( ^ = -'i2irC F M 2 B ^j2N~ c d[x] b 1 db 1 b 2 db2^ B {xiM) 

Jo Jo 

x W . [$^x 2 )<I>^x 3 ) + <J>;(x 2 )$^,(x 3 )] E<$'{tf)hf{x 1 ,X2,x 3 MM) , (A33) 

__ _ /> 1 /> OC 

M { ^ = -MttC f mI^/2N~ c d\x] hdbib 2 db 2 $ B (xi,h) 

Jo Jo 

xr/if [$^(x 2 )$^,(x 3 ) +*J(z 2 )$£-.(z 3 )] E < $\tf)hf{x 1 ,X2,x 3 M,b2) , (A34) 

M«5 = l&rC F M2,y/W c d[x] b 1 db 1 b 2 db2<S> B (x 1 ,b 1 ) 

Jo Jo 

x { [r*. (1 - x 3 )$ (x 2 )($^ (x 3 ) - (x 3 )) - r>x 2 (^(x 2 ) + ^(x 2 ))$ K . (x 3 )] 
+ [-r>(2 - x 2 )($;(x 2 ) + %(X 2 ))<S>K> (x 3 ) 

+r K ,(l+x 3 )^(x 2 )(^ K ,(x 3 )-^ s K ,(x 3 ))] Eif(tf ) )hf ) (x 1 ,x 2 ,x 3 ,b u b 2 )} , (A35) 

M% ] ab = l^C F M 2 B yj2N~ c d[x] hdhhdb^sixuh) 

Jo Jo 

x { [r x 2 ($;(x 2 ) + *%{.x 2 ))$ T K . (x 3 ) r K » (1 - x 3 )$l(x 2 )(* v K . (x 3 ) 3> a K . (x 3 ))} 
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xE { £ l {tf)h^\x u x 2l x 3l b 1 M) 

+ [r>(2 - x 2 )(<t>l(x 2 ) + %(x 2 ))<$> T K , (x 3 ) - r K , (1 + x 3 )^(x 2 )(^ K , (x 3 ) - $ a K , (x 3 ))] 
x E { a f (tf ) )hf\x 1 ,x 2 ,x 3 ,b 1 ,b 2 )} , 



Na5 



The expressions of the nonfactorizable amplitudes M. Ba and M.HeA are the same as M^ a3 and -M^ 3 but 
evolution factors E^ 3 and £^3?' replaced by E$' and respectively. 
The evolution factors are given by 

£&>'(*) =a.(t)aP'(t)S(t)\ ba=bl , 
E [ S'{t) =a s (t)a^'(t)S(t)\ b3=b2 , 

with the Sudakov factor S = S b S<pSk»- The Wilson coefficients a appearing in the above formulas are 



a 



a 



a 



a. 



a x = 

(<?)' _ 

3 — 

(<?)' _ 

4 — 

(<?)' _ 

5 — 

(<?)' _ 



Ci 
N c 
1 

n c 
i 

N~c 
1 

1 

K. 



— \C 3 + 7;e q C 9 



The hard functions 



1 and 2, are written as 



with the variables, 



,0) 



[0(h - b 2 )K Q {DM B h) I (DM B b 2 ) 

+6(b 2 - h)Ko (DM B b 2 ) I (£>M B &i)] 
xKo(DjM B b2) , forD?>0 



:|tf«(^|M s & 2 ) , forDf <0, 



h / - 2 



fl(6i - fc)^ (FMsh) Jo (FM B b 2 ) 

+6{b 2 - b!)H^ (FM B b 2 ) Jo {FM B b\) 
xKoiFjMeh) , for F? > , 



for F? < 



D 2 = 
D\ = 

Dl = 
F 2 = 

Ft = 
F 2 = 



X1X3 , 

(xi - x 2 )x 3 , 

-(1 - xi- x 2 )x 3 , 

x 2 (1 - x 3 ) , 

(xi - x 2 ) (1 - a; 3 ) , 

x\ + x 2 + (1 - xi - x 2 ) (1 - x 3 ) 



The hard scales are chosen as 



12 



tf = max ( DM B , J\D%\M b , l/h, l/h ) , 



4 X) = max [DMb, J\D\ \M b , l/h, l/h) , 



4 1] = max I FMb, \l \Fj\Mb, l/h, l/h I , 



t ( , 2) = max FM B , J \F%\M B , l/h, l/h ) • (A43) 
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TABLE I. Helicity amplitudes and relative phases. 



Mode 


BR(W- a ) 




kin r 




<j>\\ (rad.) 


4>±(rad.) 


07T U 


14.86 


0.750 


0.135 


0.115 


2.55 


2.54 



JK^ 15^96 0748 0133 am 2l35 2^54 



TABLE II. Helicity amplitudes and relative phases: (I) without annihilation and nonfactorizable contributions, (II) without 
annihilation contributions, and (III) without nonfactorizable contributions. 



Mode 


BR{lQ- b ) 






\A±\" 


011 (rad.) 


4>±(rad.) 


0if* u (I) 


14.48 


0.923 


0.040 


0.035 


IT 


TV 


(II) 


13.25 


0.860 


0.072 


0.063 


3.30 


3.33 


(III) 


16.80 


0.833 


0.089 


0.078 


2.37 


2.34 


<j>K*+(l) 


15.45 


0.923 


0.040 


0.035 


7T 


7T 


(II) 


14.17 


0.860 


0.072 


0.063 


3.30 


3.33 


(III) 


17.98 


0.830 


0.094 


0.075 


2.37 


2.34 
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